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Abstract 



We give a general analysis of neutrino mixing in the seesaw mechanism 
with three flavors. Assuming that the Dirac and u-quark mass matrices are 
similar, we establish simple relations between the neutrino parameters and 
individual Majorana masses. They are shown to depend rather strongly on the 
physical neutrino mixing angles. We calculate explicitly the implied Majorana 
mass hierarchies for parameter sets corresponding to different solutions to the 
solar neutrino problem. 
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I. INTRODUCTION 



One of the most pressing questions in particle physics has been the determination of 
the intrinsic properties of the neutrinos, namely, their masses and mixing angles. The 
recent atmospheric neutrino data |L| suggest strongly that neutrinos do have masses, and 
that, unlike the quark sector, at least some of the mixing angles are large. The most 
appealing model for small neutrino masses derives from the seesaw mechanism ||. In doing 
so, however, one also introduces additional unknowns in the form of the Dirac and Majorana 
mass matrices. We do have a handle on the Dirac matrix, since, from the ideas of GUTs, it 
should be similar to that of the quark sector. Not much is known about the Majorana mass 
matrix. The challenge is then to find out what the Majorana matrix is like in order for the 
effective neutrino mass matrix to come out correctly. 

In this paper, we will analyze the general structure of the seesaw mass matrix. Without 
loss of generality, we will work in the basis in which the charged lepton and Majorana mass 
matrices are diagonal. As a starting point, we assume that the Dirac mass matrix, in analogy 
to the quark mass matrix, has hierarchical eigenvalues and small left-handed mixing angles. 
Even in this case, large mixing can occur through the interplay of the Dirac and Majorana 
matrices 0-0. 

In the seesaw mechanism, the effective neutrino mass matrix is given by 



m cS = m D M L m 



D 



(1) 



The Dirac matrix can be written as 



m D = U m% as V 



(2) 



in the basis where M 1 is diagonal, 
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Here Uq and Vq are left-handed (LH) and right-handed (RH) rotations, respectively. m^ ag 
a diagonal matrix with eigenvalues m« (i = 1, 2, 3). The Majorana mass eigenvalues are given 
by Mi — 1/Ri (i — 1, 2, 3). For simplicity, we assume that Uq and Vq are real (i.e. we ignore 
CP violating effects). Note that Vq also contains contribution from the diagonalization of 
M -1 , and thus it may contain large angles. However, we will restrict our discussions to 
small angles in Vo- 
lt is convenient to write 



U^w^Uq = NN 1 



(4) 



where 



N = m^VoM- 1 ! 
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For hierarchical Dirac masses, to a good approximation 
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Here = (Vb)y are the matrix elements of Vb- More precisely M, the leading cor- 
rection to this approximation is a LH rotation with rotation angles (0125 0135 023) — 
/rrnlii a^a alk] _ This can be absorbed into Un, and will be ignored henceforth. 

V mi V22 ' "13 V33 ' raz V33 / ' 

If M~ l oc I, then NN T will be diagonal and m eff will be diagonalized by Uq. When M -1 
deviates from the identity matrix, the product VqM~2 is no longer orthogonal, then 



N = UN di ^W 



(7) 



where U and W are LH and RH rotations respectively. The LH rotation U is induced by 
so that 



^Uq 1 ) m eff (U U) = (A^ diag ) 2 



(8) 



Thus, the mixing angles for the effective neutrino mass matrix come from UqII, while the 

eigenvalues of iV are just ^m| ff , with mf 1 denoting the effective neutrino masses. To the 
extent that we may assume Uq to be nearly the identity matrix, the case of large mixing 
angles receives its main contribution from U, which is induced by M. 



From Eq. (pH, we have 
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(9) 



While we have assumed that the Dirac mass eigenvalues m; have the hierarchical structure 
of the quark masses, little is known about the values of Ri. 

The problem within the seesaw model is then twofold. One, given the parameters m;, 
Ri, and Vy, what is the LH rotation matrix U as defined in Eq. (0) ? Conversely, if we 
know U, i.e. the physical mixing angles for the neutrinos, what can we deduce about these 
parameters? We will discuss these issues in the following sections. 



II. GENERAL ANALYSIS 

Let us start from a general 3x3 matrix 
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Anticipating the applications to the neutrino sector, we will assume 

a 2 < b 2 ~ c 2 , (11) 

where we have used the notation a = (ai,a 2 ,a 3 ), etc. We have shown elsewhere [|J that 
when a matrix is brought into the upper triangular form, we can obtain the LH mixing angles 
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easily. Now, it is always possible to find a RH rotation so that N becomes upper triangular. 
In fact, geometrically, this amounts to a new coordinate system where c is aligned with 
the third axis, while the second and first axes are in the directions c x (b X c) and b x c, 
respectively. It is then clear that, for an appropriate RH rotation R, 




N t = NR= \b x c\ b.c , (12) 



where i = bx c/\bx cj, j = cx (b x c) / (\c\\b x c\) , and c = c/\c\. Note that, since the physical 
mass matrix is given by NN T , the matrix N in Eq. (|10|) is arbitrary up to any RH rotation. 
The matrix in Eq. (0), however, is constructed entirely from rotational invariants, and is 
given in terms of physical quantities. There is still some ambiguity in Eq. ([12]) , in that the 
diagonal elements can be of either sign, corresponding to the choice of orientation of the 
axes. 

To diagonalize N, we first diagonalize its (23) submatrix explicitly by Rl{23) and Rr(23), 
with 



tan 26* 



L 



2\b.c\ 



23 ~ c 2 - 6 2 

tan 20* = 2 fc^ X j . (13) 



c 4 + (b.c) 2 -(bxc) 2 
Further, the eigenvalues of N are fi 2 and /x 3 , with 



A 2 = ^(c 2 + b 2 ) 2 -\bxc\ 2 . (14) 



Thus we find 

Ni = R L (23)N t R R (23) 



f a.i «2 «3 

fi2 | , (15) 
V 0/i 3 . 



where (02,^3) = (a.j,a.c)Rji(23). Note that the mixing angle 6^ is maximal when b 2 = c 2 . 
Also, a mass hierarchy, /Z3 3> //2, implies that b.c ~ c 2 , i.e. b and c are nearly parallel. 
Since |a 3 | < a // 3 , a LH (13) rotation Rl(13), with ~ \as\/fi3 1, removes the (13) 
element of N' t without changing the other elements appreciably. Summarizing, we see that, 
if a 2 < c 2 , 

a.i ai \ 

iV = ^(13)^(23)^^(23) ~ [ fi 2 , (16) 

fi 3 J 



where Rl(13) ~ J, i? is defined in Eq. (p!2[) , and the rotation angles in R L R {23) are given 
in Eq. (|i~3"P . The final diagonalization of iV can be achieved by a combined LH and RH (12) 
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rotation. In particular, maximal mixing is possible if «2 ~ fi2- We emphasize that as long 
as the (13) rotation i?^(13) is small, the diagonalization of N can be decomposed into that 
of two 2x2 matrices. This is a very useful simplification. 

Of particular interest in neutrino physics is the possibility of large mixing angles. As- 
suming \a\ <C |c|, which ensures that 0f 3 <C 1, maximal (23) rotation, 9% 3 = ir/4, is obtained 
if \b\ = |c|. There are now two possibilities: 

A) \b.c\ > \b x c] 2 . 

In this case, from Eq. (|TJj) we find that /z 2 ~ \b x c| 2 /(c 2 + b 2 ) <C /i 2 . Under this condition, 
we can have bimaximal mixing, 9\ 2 ~ 7r/4. This is achieved if a 2 — [i2- Note that, for 
|6 x c| \b-c\, we have the bound |tan26 |2 ? 3 | < \b x c|/|6.c|. Thus | sin^ffa.cl <C /i2, so the 
condition a 2 — /i2 can also be written as a.j ~ /i 2 . 

B) \b.c\ ~ |6 x 5|. 

In this case /12 ~ ^3. All of the components of a are smaller than i n particular, | /X2 1 ^> &2- 
Thus, 9 12 ~ \oi2 |/|a*2 I "C 1, and we can only have single maximal mixing. 

Knowing N, the above analysis yields the neutrino mixing angles Ok. In Eq. (|S|), N 
is given in terms of the parameters rrii, Ri, and V^. There are generally two classes of 
possibilities: (A) All Ri are of the same order, or (B) there is a strong hierarchy in the 
Majorana sector. 

For the case of no Majorana hierarchy, R\ ~ R2 ~ R3, we can find approximately the 
rotation angles, which give 

fc-MS^Ws^MS*)- (17) 

so that 

R L NN T R T L = iV 2 ag . (18) 
In other words, the induced neutrino mixing angles are given approximately by 

Of j * ^f/ji > J) = (!> 2), (2, 3), (1, 3). (19) 

Unless there is a stronger hierarchy in the Majorana sector in comparison to the Dirac sector, 
all neutrino mixing angles are negligible when Vij is reasonably small. 

To analyze the situation when there is a strong hierarchy in Ri, we will limit our discussion 
to situations which are physically interesting. For this purpose, it is actually more convenient 
to reverse our procedure above, i.e., given the physical mixing angles, we deduce the form of 
the matrix N. We note that, experimentally, it is quite suggestive that the neutrino masses 
are hierarchical, and that the neutrino mixing matrix U a i, (a = e,/i, r;i = 1,2,3), is given 
approximately by: 

U = R 23 ((t>)Ri3(e)Ri2(9) 

( c e s e 

— -{s$C4, + ec e s^) (cec^ - eses^) . (20) 
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where tan ~ (9(1) and e Cl 0. We have kept only the leading terms of e in U. Here, tan 9 
can be either of order 1 (corresponding to the "bi- maximal" mixing scenario), or could be 
smaller ("single-maximal" scenario). In the latter situation, we will discuss two possibilities: 
(A) e <C sq, (B) e ^> sq. Appropriate approximations can be obtained for each of them. 
For the case eCsj, the matrix N is 





( n% 







N = U 





n 2 


•)-( 
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c e ni s e n 2 en 3 
-secern c e c^n 2 s^n 3 
ses^ni -c e s^n 2 c^rij, 



up to an arbitrary RH rotation. ^From Eq. 
readily seen that 



(21) 



in, with the assumption ni < n 2 <C «3, it is 



c2 



= Tin 



m 



eff 



(22) 



and 



\b x c| c e n 2 n 3 c e n 2 

IM 



< 1 . 



(23) 



These two equations imply that b and c are nearly parallel to leading order in n 2 /n 3 . Also 
a <C b,c, so that our general analysis is applicable to N. 

We can put Eq. (|21|) in the lower triangular form through a RH rotation. We have 



N 






(axb).c 
\axb\ J 



fa 

b. a/a \a x b\/a 

c. a/a 

\ ' a|axb| 

/ yfifyfi + s 2 e n 2 2 + e 2 nl 

(c ( j } sec e nl + s^enD/a ses^n^/a • (24) 

\(-s 4l s e c e nl + c^enfj/a s e c^n 2 n 3 /a nx/sgs^J 

is uniquely determined by the physical 







Note that, as in Eq. (|T2"D, the entries of N in Eq. 
neutrino parameters, while Eq. fl2~l"|) is not. 

A similar analysis can be done for e ^> sg. Here, Eqs. 
To lowest order in sg and e and with = 7r/4, we obtain 



lj) and (|24 ) are no longer valid. 



N 



I tn 3 
vl( 1 + 












(25) 



We may now compare Eq. (^j) to Eqs. (|24|) and (|25|). Given the neutrino parameters, n*, 
Sg, and e, this method gives immediately the Majorana masses Mi, assuming that we may 
identify the Dirac masses with the quark masses. We emphasize that this comparison is 
viable because the lower triangular form for the N matrix is unique, there being no more 
ambiguities due to RH rotations. Also, the structure in these equations shows immediately 
that they are consistent only if there is a very strong hierarchy, Ri 3> R 2 3> -R3. In the next 
section, we will discuss in detail the physical consequences of this result. 
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III. DISCUSSIONS AND NUMERICAL RESULTS 



Assuming hierarchical Dirac neutrino masses and small RH rotation angles, it is seen 
from Eqs. (0) that the N matrix is naturally of lower-triangular form. The discussion of 
physical constraints can be done most conveniently in this lower triangular basis. In this 
section, we will start from Eqs. @, (0), and ( pop and present the results for the neutrino 
parameters derived from different solutions to the solar and atmospheric neutrino observa- 
tions. To be concrete, we will use = 45°, and = (m u ,m c ,m t ) (at the seesaw scale), 
although different choices can be accommodated. 



(A) e < s e : 



Depending on e, Eq. (p4|) gives rise to three different patterns, with the largest elements 
located in the (a) (2,2), (3,2), (b) i = 1,2, (c) (2,1), (3,1) positions, respectively. It 

can be seen that only type (a) with e <C mf /ui 3 is natural, which we will concentrate on. 
In this limit, with tan# > rii/n 2 , 



( s e n 2 



N 



i 



V2 
1 

n/2 



cen 2 




(26) 



where we identify n 2 



m 



eff 



Comparison to Eq. (ffi) immediately yields (with Vn ~ Vj 



22 



V, 



33 
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eff 



m: 



rrin 



,cff 



m 



eff 



n-. 



I = 2(R 2 m 2 ) 2 
s 2 e (R 3 m 3 ) 2 /2 



sjM 1 
2m 2 c 

2M, 



(27) 



_]^ = V^ = m 2 = mc 
V 21 V 22 m 3 m t ' 

Quite surprisingly, m| scales as m 2 rather than m 2 . This is because in Eq. (|26|) the (22) 
element is one of the largest. This gives a scale for M 2 much lower than one would expect, 

2m 2 

;6xl0 9 GeV, (29) 



Am 



2 



where we have used m c (M 2 ) ~ 0.4 GeV, and mf 1 = yjAm 2 tm = V3 x 10~ 3 eV, which will 
also be used in the following. 

Note also that V 32 is independent of the Majorana mass Mi and that it scales linearly with 
m c /m t . This means that the Majorana sector decouples owing to its very large hierarchy. 
Similarly, we have 
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^= r mi = r mu , (30) 
hi v 2 tan 61712 v 2 tan #m c 

V31 mi -m u 



V11 a/2 tan #m 3 a/2 tan 



(31) 



^From Eqs. (H), (i^)' (0)) we see that a ^ °f the angles are small. 

The other two heavy Majorana mass values depend on the different solutions to the 



solar neutrino problem [ 10| . We will use m| ff = yAm^m, wijp = yj Am s 2 olar , while mf 1 is 
not known. However, we can derive a bound for M3 using the parameter r = mf /mf ^> 1. 
From Eq. (p7[), we note that M x and M 3 , besides the usual Dirac mass squared, depend 
sensitively on sq as well as on the effective neutrino masses. This can be seen directly in the 
following equations: 

Mi m\mf M 2 _ Am 2 c mf 

~W 2 ~ 2m 2 c s 2 e mf ' M3 ~ m 2 t s 2 e mf ' ^ ' 

We will now turn to numerical estimates with inputs coming from the three solutions to 
the solar neutrino problem, vacuum oscillations (VO), large angle MSW (LAM), and small 
angle MSW (SAM). 

Al. VO 



The neutrino mass and mixing can be taken as [I0|, 



6 ~ 45°, mf = \JAM 2 olav = V7 x 10" 11 eV. (33) 
^From Eq. (|32|), we find 

2m 2 

Mi = — ^ « 5 x 10 8 GeV , (34) 
m 2 

M 3 /r « 4 x 10 17 GeV (r = mf /mf > 1) . (35) 



A2. LAM 
Here we take ITO 



sin 2 2# solar ~ 0.8 , mf « y Am s 2 olar = V3 x 10" 5 eV . (36) 
Going through the same analysis as in the VO case, we have: 

Mi«lx 10 6 GeV , M 2 « 6 x 10 9 GeV , M 3 /r « 4 x 10 14 GeV (r > 1) . (37) 



S 



A3. SAM-I 



We now have [ID 



sin 2 29, 



solar 



5 x 10 



-3 



i n':," ks \/Am 2 olar = V5 x 10~ 6 eV 



(38) 



They imply 



M!«7x 10 8 GeV , M 2 « 6 x 10 9 GeV , M 3 /r « 4 x 10 12 GeV (r > 1) . (39) 

The above results were obtained under the assumption that e — > 0. Thus, if e ^> ^2/^3, 
Eq. (|26|) is no longer valid. Also, if e ^> sg, Eq. (|24]) should be replaced by Eq. (^). This 
condition is actually very likely to be valid for the case of the SAM solution. We treat this 
case in detail next. 



(B) e > s e (SAM-II): 

The parameters for 9 and m| ff are taken as in SAM-I. However, we now use Eq. fl25|) instead 
of Eq. (f?|). Also, the numerical result depends on the value of e, which, for defmiteness, we 
will take to be e = 0.1. We obtain 



M^2x 10 6 GeV , M 2 «lx 10 11 GeV , M 3 /r « 2 x 10 13 GeV (r > 1) . (40) 



In summary, we find that if the physical neutrino parameters are known, we can obtain 
the Majorana masses directly when we identify the Dirac masses with the quark masses. 
The Majorana masses have rather strong dependence on the physical mixing angles, so that 
their values span a wide range. Numerically, it is noteworthy that M 2 ~ 6 x 10 9 GeV for a 
wide range of parameters. Also, the VO solution for M 3 (3> 4 x 10 17 GeV) seems too large 
for it to be viable. Finally, in some cases Mi can be rather low (~ 10 6 GeV). 



IV. CONCLUSION 



The observation that neutrino masses are tiny has a natural explanation in the seesaw 
model. However, its complicated structure also means that the neutrino mixing angles do 
not derive simply from the seesaw components, viz., the Dirac matrix mp and the Majorana 
matrix M~ l . In fact, if we write m D = f/ m* ag F , M -1 = Um{M~^/^) 2 U'^^ m e s depends, 
roughly speaking, quadratically on all of the components that we displayed. 

In this paper, we analyze the problem in several steps. First, we write m e ^ = NN T , 
so that iV depends linearly on the aforementioned components. But iV has the further 
ambiguity of an arbitrary RH rotation. We will eliminate this ambiguity by reducing iV to the 
triangular form. The lower triangular form arises naturally if vtld has a hierarchy. However, 
the upper triangular form is the easiest to use in order to extract the LH, physical, neutrino 
mixing angles. By expressing the triangular matrix elements in terms of RH rotational 
invariants, it is easy to transform back and forth between the different forms of N. Thus, 
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given the parameters in mo and M, we can deduce the neutrino mixing angles. Conversely, 
given the physically plausible values of the neutrino masses and mixing angles, we can obtain 
the constraints that must be satisfied by the parameters in and M. Experimentally, 
we have a fairly good idea about the intrinsic neutrino parameters. Their masses are most 
probably hierarchical, the (23) mixing angle is almost maximal, and the (13) mixing angle 
is small. From these parameters, we can infer the properties of the Majorana masses and 
the RH mixing angles of tud, if we assume that is similar to the u-quark mass matrix, 
i.e., rrt£, has small LH mixing angles and mj ~ (m u ,m c ,m t ). It was found that there must 
be a large hierarchy in M, proportional not only to the ratios of Dirac masses squared, but 
also to the squares of the mixing angles. In addition, the RH mixing angles in combined 
with M _1 must be very small, and must be equal to the mass ratios in mo- Physically, the 
large hierarchy in M implies the existence of intermediate mass scales. It would be most 
interesting if these conclusions can be corroborated by other sources. 

Using solar neutrino solutions as inputs, we calculated the individual Majorana masses. 
Because of their strong dependence on the physical neutrino mixing angles, a wide range of 
values was found. In particular, M3 is so large (^> 4 x 10 17 GeV) for the VO solution which 
makes it highly disfavored. In this work we have not treated the issue of renormalization, 
although it can be shown that the RGE effects are small [O. We have also not discussed 



the case when either or both of mo and M are complex. Although the 2x2 problem can 
be solved, the 3x3 case does not seem to have a simple solution. Nevertheless, it can be 
shown that, with hierarchical masses and small angles, complex phases do not contribute 
significantly 0]. We hope to return to this problem in the future. 
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